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The fast diffusion equation u t = (u u x )x is investigated from the symmetry point of view in 
development of the paper by Gandarias [Phys. Lett. A 286 (2001) 153-160]. After studying equiv- 
alence of nonclassical symmetries with respect to a transformation group, we completely classify 
the nonclassical symmetries of the corresponding potential equation. As a result, new wide classes 
of potential nonclassical symmetries of the fast diffusion equation are obtained. The set of known 
exact non-Lie solutions are supplemented with the similar ones. It is shown that all known non-Lie 
solutions of the fast diffusion equation are exhausted by ones which can be constructed in a regular 
way with the above potential nonclassical symmetries. Connection between classes of nonclassical 
and potential nonclassical symmetries of the fast diffusion equation is found. 



1 Introduction 

Investigation of nonlinear heat (or diffusion if u represents mass concentration) equations by 
means of symmetry methods was started as early as in 1959 with Ovsiannikov's work [19] where 
the author performed the group classification of the class of equations of the form 



Nonclassical symmetries of equations from class (pQ) were investigated in [3, 10, 11, 15]. In 
particular, the authors of [11] obtained the determining equations for the coefficients of condi- 
tional symmetry operators for the wider class of nonlinear reaction-diffusion equations of the 
form ut = (f(u)u x ) x + g{ u ) and constructed a number of their exact solutions. Review of results 
on symmetries, exact solutions and conservation laws of such equations is given e.g. in [16]. 

The diffusion processes described by (pQ) are known to arise in different fields of physics 
such as plasma physics, kinetic theory of gases, solid state and transport in porous medium. 
In many metals and ceramic materials the diffusion coefficient f{u) can, over a wide range of 
temperatures, be approximated as u~ a , where < a < 2 [29]. So, one of the mathematical 
model of the diffusion processes is 



Equations ([2]) are called fast diffusion equations in the case < a < 2 since these values of a 
correspond to a much faster spread of mass than in the linear case (a = 0). 

In this Letter we restrict ourselves with the special case a = 1, i.e. with the equation 





ut = (u a u x ) x . 



(2) 



u t = (u u x ) x . 



(3) 



It emerges in plasma physics as a model of the cross-field convective diffusion of plasma including 
mirror effects and in the central limit approximation to Calerman's model of the Boltzmann 



equation. Equation ([3]) governs the expansion of a thermalized electron cloud described by 
isothermal Maxwell distribution. It is also the one-dimensional Ricci flow equation. (See [4,5,29] 
and references therein.) 

Equation ([3]) has a number of remarkable mathematical properties which distinguish it from 
class ([2]). Thus, ([3]) can be rewritten in the form ut = (lnu) xx whereas for the other values of a 
the function under d xx is a power one. It admits a discrete potential invariance transformation. 
For this equation wide classes of exact solutions were constructed in a closed form while reduction 
of ([5]) in the general case results in ordinary differential equations which usually cannot be 
integrated explicitly. Its potential form admits two kinds of variable separation. 

The fact that ([3]) is written in a conserved form allows us, following Bluman et al. [7-9], to 
consider the corresponding potential system 

v x = u, v t = u~ x u x (4) 

and to find potential symmetries of equation (|3|). Namely, any local symmetry of system Q 
induces a symmetry of the initial equation ([3]) . If transformations of some of the "non-potential" 
variables t, x and u explicitly depend on the potential v, this symmetry is a nonlocal (potential) 
symmetry of equation ([3]). 

It follows from that the potential v satisfies the nonlinear filtration equation 

Vt = V^Vxx (5) 

with the special value v x ~ l of the filtration coefficient. We will also call equation ([5]) the potential 
fast diffusion equation. Akhatov, Gazizov and Ibragimov carried out group classification of the 
nonlinear filtration equations of the general form 

vt = f(v x )v xx (6) 

and investigated their contact and quasi-local symmetries [1,2,16]. 

Lie symmetries of ([3]) are well known (see Section 2). All its exact solutions constructed in 
closed form by reduction with Lie symmetries are listed e.g. in [21]. 

Some non-Lie exact solutions of ([3]) were obtained in [15,28,29]. Thus, Rosenau [29] found 
that equation ([5|) admits, in addition to the usual variable separation v = T(t)X(x), the additive 
one v = Y(x + Xt) + Z(x — At) which is a potential additive variable separation for equation ([3]). 
To construct nonclassical solutions of © , Qu [28] made use of generalized conditional symmetry 
method, looking for the conditional symmetry operators in the special form Q = (u xx +H(u)u x 2 + 
F(u)u x + G{u))d u . Gandarias [15] investigated some families of usual and potential nonclassical 
symmetries of ([2]). In particular, using an ansatz for the coefficient rj, she found nontrivial 
reduction operators in the so-called "no-go" case [12,22,30] when the coefficient of dt vanishes, 
i.e. operators can be reduced to the form Q = d x + r/(t, x, u)d u . 

In the recent paper [9] a preliminary analysis of nonclassical symmetries of equations 
from class ([6]) was performed. A more detailed consideration was carried out for the case 
/ = {v 2 x + v x ) , and only some examples of reduction operators and corresponding exact so- 
lutions were constructed. Let us note that equation ([6]) with / = (v x + v x )~ l is reduced by 
the point transformation t = t, x = x + v, v = v to equation © which corresponds to the 
value / = v x ~ x and is simpler and more convenient for investigation. All results on symme- 
tries and exact solutions of the equation from [9] can be derived from the analogous results for 
equation ([5]). 

In this Letter the fast diffusion equation ([3]) is investigated from the symmetry point of view. 
The nonclassical symmetries of the corresponding potential equation ([5]) are completely classified 
with respect to its Lie symmetry group. As a result, new wide classes of potential nonclassical 
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symmetries of equation ([3]) are found. Some classes of potential nonclassical symmetries prove 
to be connected with usual nonclassical ones on the solution set of potential system (j3J). The 
set of exact non-Lie solutions constructed in [15, 28, 29] is supplemented with the similar ones. 
It is shown that all known non-Lie solutions of the fast diffusion equation are exhausted by ones 
which can be constructed with the above potential nonclassical symmetries. 

Our Letter is organized as follows. First of all (Section [2]) we adduce results on Lie and 
potential symmetries of ([3|), including discrete ones. It is important since classical symmetries 
really are partial cases of nonclassical symmetries and below we solve the problem on finding 
only pure nonclassical symmetries which are not equivalent to classical ones. Moreover, our 
approach is based on application of the notion of equivalence of nonclassical symmetries with 
respect to a transformation group, which is developed and investigated in Section O Usage of 
equivalence with respect to the complete Lie invariance group including the discrete symmetries 
plays a significant role in simplification of proof, testing and improving presentation of the main 
result (Theorem [TJ Section H]). In spite of the techniques applied in [15] and similarly to [9], 
we use the single potential equation (0) instead of potential system @, to produce potential 
nonclassical symmetries of ([3]). After the "no-go" case of the zero coefficient of dt is discussed, 
all the reduction operators having the nonvanishing coefficient of dt are classified. Connection 
between partial classes of usual and potential reduction operators of (pQ) is studied in Section OJ 
In Section [6] the known Lie solutions of ([3]) and ([5]) are collected. A list of non-Lie solutions 
is supplemented with the similar ones. Connections between exact solutions and different ways 
of their construction are discussed shortly. In conclusion some recent results on nonclassical 
symmetries of equations (J6|) are announced. 

2 Lie and potential symmetries of fast diffusion equation 

The Lie invariance algebra 

Ai = (d t , d x , td t + ud u , xd x - 2ud u ) 

of equation ^ was found in [19]. The complete Lie invariance group G\ of ([3J) is generated 
by both continuous one-parameter transformation groups with infinitesimal operators from A\ 
and two involution transformations of alternating sign in the sets {t, u} and {x}. Action of any 
element from G\ on the function u is given by the formula 

u(t, x) = £4 u{e^t + e\, e 4 x + e 2 ), 

where e\, . . . , £4 are arbitrary constants, £364 7^ [19]. 

The Lie symmetry properties of ([3]) are common for diffusion equations. Uncommonness of 
equation ([3]) from the symmetry point of view becomes apparent after introducing the poten- 
tial v and considering potential system ^ or potential equation ([5]). Point and nonclassical 
symmetries of flD) or ([5]) are called potential and potential nonclassical symmetries of © corre- 
spondingly. 

The Lie invariance algebra 

A 2 = (d t , d x , d v , td t + vd v , xd x - vd v ) 

of equation ^ and the corresponding connected Lie symmetry group are quite ordinary for 
nonlinear filtration equations. However, equation ([5]) is distinguished for its discrete symme- 
tries since it possesses, besides two usual sign changes in the variable sets {t,v} and {x,v}, 
the hodograph transformation t = t, x = v, v = x. These three involutive transformations to- 
gether with the continuous one-parameter transformation groups having infinitesimal operators 
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from A2 generate the complete Lie invariance group G2 of ([5]). Therefore, G2 consists of the 
transformations 

t = S3t + si, x = e^x + £2, v = £3£4 -1 t> and 

t = £^t + ei, X = E3E4~ 1 V, V = E4X + £2, 

where £1, . . . , £4 are arbitrary constants, £3£4 7^ 0. 

A similar result is true for system Q. Namely, it is invariant with respect to the following 
transformation 

t = t, x = v, u = u~ 1 , v = x (7) 

which is additional to the usual Lie symmetry group G\ of equation © and is called the potential 
hodograph transformation of this equation. 

It can be proved [25] that the set of Lie invariant solutions of equation ([3j) is closed under 
transformation ([7|). 

3 Equivalence of reduction operators 
with respect to transformation groups 

The notion of nonclassical symmetry was introduced in 1969 [6]. A precise and rigorous definition 
was suggested later (see e.g. [13,31]). 

Consider an rth order differential equation C of the form L(t,x,ur r \) = for the unknown 
function u of two independent variables t and x. Here ut r \ denotes the set of all the derivatives 
of the function u with respect to t and x of order not greater than r, including u as the derivative 
of order zero. Within the local approach the equation C is treated as an algebraic equation in 
the jet space of order r and is identified with the manifold of its solutions in J^: 

C = {(t, x, U( r )) € Jv> I L(t, x, ii( r )) = 0}. 

The set of (first-order) differential operators of the general form 

Q = r(i, x, u)d t + £(*, x, u)d x + rj{t, x, u)d u , (r, £) / (0, 0), 

will be denoted by Q. Here and below dt = d/dt, d x = d/dx and d u = d/du. Subscripts of 
functions denote differentiation with respect to the corresponding variables. 

Two differential operators Q and Q are called equivalent if they differ by a multiplier being 
a non-vanishing function of t, x and u: Q = XQ, where A = X(t, x, u), X ^ 0. The equivalence of 
operators will be denoted by Q ~ Q. Denote also the result of factorization of Q with respect 
to this equivalence relation by Qf. Elements of Qf will be identified with their representatives 
in Q. 

The first-order differential function Q[u] := r](t,x,u) — r(t,x,u)ut — £(t,x,u)u x is called the 
characteristic of the operator Q. The characteristic PDE Q[u] =0 (called also the invariant 
surface condition) has two functionally independent integrals x, u) and u(t, x, u). Therefore, 
the general solution of this equation can be implicitly presented in the form F(£, uj) = 0, where F 
is an arbitrary function of its arguments. 

The characteristic equations of equivalent operators have the same set of solutions. And vice 
versa, any family of two functionally independent functions of t, x and u is a complete set of 
integrals of the characteristic equation of a differential operator. Therefore, there exists a one-to- 
one correspondence between Qf and the set of families of two functionally independent functions 
of t, x and u, which is factorized with respect to the corresponding equivalence. (Two families 
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of the same number of functionally independent functions of the same arguments are considered 
equivalent if any function from one of the families is functionally dependent on functions from 
the other family.) 

Since (r, £) ^ (0,0) we can assume without loss of generality that Q u ^ and ^ and 
resolve the equation F = with respect to £: Q = ip(u>). This implicit representation of the 
function u is called an ansatz corresponding to the operator Q. 

Denote the manifold defined by the set of all the differential consequences of the characteristic 
equation Q[u] =0 in by Q^ r \ i.e. 

Q (r) = {(i,x,u (r) ) € J (r) \D?D%Q[u] = 0, a,/3 G NU{0}, a + (3 <r}, 

where D t = dt + u a+ \^d Uaf3 and D x = d x + u a ^ + id Ual3 are the operators of total differentiation 
with respect to the variables t and x, the variable u a p of the jet space corresponds to the 
derivative d a+ ^u/dt a dx^ . 

Definition 1. The differential equation C is called conditionally invariant with respect to 
the operator Q if the relation Q^L(t,x, u (r))\£ ri Q(r) = holds, which is called the conditional 
invariance criterion. Then Q is called an operator of conditional symmetry (or Q-conditional 
symmetry, nonclassical symmetry, etc.) of the equation C. 

In Definition [T] the symbol Qr r \ stands for the standard rth prolongation of the operator Q 

[18,20]: Q (r) =Q + Eo<a+/3<r ^ d u a0 , where rfP = D?D P x Q[u] + tu q+1j/3 + £« a ,/?+i- 

The equation C is conditionally invariant with respect to the operator Q iff the ansatz con- 
structed with this operator reduces C to an ordinary differential equation [31]. So, we will also 
call operators of conditional symmetry by reduction operators of C. 

Lemma 1 ( [14,31]). If the equation C is conditionally invariant with respect to the operator Q, 
then it is conditionally invariant with respect to any operator which is equivalent to Q. 

The set of reduction operators of the equation C is a subset of Q and so will be denoted by 
Q(C). In view of Lemma [TJ Q € Q(£) and Q ~ Q imply Q G Q(£), i.e. Q(C) is closed under 
the equivalence relation on Q. Therefore, factorization of Q with respect to this equivalence 
relation can be naturally restricted on Q(C) that results in the subset Qf(C) of Qf. As in the 
whole set Qf, we identify elements of Qf (£) with their representatives in Q(C). In this approach 
the problem of complete description of reduction operators for the equation C is nothing but 
the problem of finding Qf(£). 

We can essentially simplify and order classification of reduction operators, additionally taking 
into account Lie symmetry transformations of an equation or equivalence transformations of a 
whole class of equations. 

Lemma 2. Any point transformation of t, x and u induces a one-to-one mapping of Q into 
itself. Namely, the transformation g: t = T(t,x,u), x = X(t,x,u), u = U(t,x,u) generates the 
mapping g*: Q — > Q such that the operator Q is mapped to the operator g*Q = fd^ + £d x + fjdy,, 
where f(t,x,u) = QT(t,x,u), £(t,x,u) = QX(t,x,u), fj(t,x,u) = QU(t,x,u). If Q' ~ Q then 
9*Q' ~ g*Q- Therefore, the corresponding factorized mapping g^: Qf — > Qf also is well-defined 
and one-to-one. 

Definition 2 ( [23,26]). The differential operators Q and Q are called equivalent with respect 
to a group G of point transformations if there exists g € G for which the operators Q and g*Q 
are equivalent. Notation: Q ~ Q mod G. 

Lemma 3. Given any point transformation g of the equation C to an equation C, g* maps Q(C) 
to Q{£) in a one-to-one manner. The same statement is true for the factorized mapping gi from 
Q f (£) to Q f (£). 
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Corollary 1. Let G be a Lie symmetry group of the equation C. Then the equivalence of operators 
with respect to the group G generates equivalence relations in Q(C) and in Q{{C). 

Consider the class C\s of equations Cg: L(t,x,ur r \,9(t,x,ur r ))) = parameterized with the 
parameter-functions 9 = 9(t,x,ur r \). Here L is a fixed function oft, x, ur r \ and 9. 9 denotes the 
tuple of arbitrary (parametric) functions 9(t,x,u^) = (9 1 (t,x,u^), . . . ,9 k (t,x,u^)) running 
the set S of solutions of the system S(t,x,u^,9^(t,x,u^)) = 0. This system consists of 
differential equations on 9, where t, x and U( r ) play the role of independent variables and 9^ 
stands for the set of all the partial derivatives of 9 of order not greater than q. In what follows 
we call the functions 9 arbitrary elements. Denote the point transformations group preserving 
the form of the equations from C\s by G~. 

Let P denote the set of the pairs each of which consists of an equation Cg from C\s and an 
operator Q from Q{Cg). In view of Lemma El action of transformations from the equivalence 
group G~ on C\s and {Q(Cg) \ 9 £ S} together with the pure equivalence relation of differential 
operators naturally generates an equivalence relation on P. 

Definition 3. Let 9,9' eS,Q G Q(£g), Q' € Q{C e >). The pairs (Cg, Q) and (C e >,Q') are called 
Gr~ -equivalent if there exists g G G~ which transforms the equation Cg to the equation Cg', and 

Q' ~ g*Q- 

Classification of reduction operators with respect to G~ will be understood as classification 
in P with respect to the above equivalence relation. This problem can be investigated in the 
way that is similar to usual group classification in classes of differential equations. Namely, we 
construct firstly the reduction operators that are defined for all values of the arbitrary elements. 
Then we classify, with respect to the equivalence group, the values of arbitrary elements for each 
of that the equation Cg admits additional reduction operators. 

In an analogues way we also can introduce equivalence relations on P, which are generated 
by either generalizations of usual equivalence groups or all admissible point transformations [24] 
(called also form-preserving ones [17]) in pairs of equations from C\s. 

4 Reduction operators 

of nonlinear filtration equation 

In this section we describe GVinequivalent reduction operators of the potential fast diffusion 
equation ©. Here reduction operators have the general form Q = rdt + £d x + 9d v , where r, £ 
and 9 are functions of t, x and v, and (r, £) ^ (0,0). Since ([5]) is an evolution equation, there 
are two principally different cases of finding Q: t ^ and r = 0. 

In the case r = we have ^ ^ 0, and up to the usual equivalence of reduction operators we 
can assume that £ = 1, i.e. Q = d x + 9d v . The conditional invariance criterion implies only one 
determining equation on the coefficient 9 

99t = 9 XX + 299 xv + 9 2 9 vv — 9 1 (9 X ) 2 — 29 X 9 V — 9{9 V ) 2 

which is reduced with a non-point transformation to equation ([5]), where 9 becomes a parameter. 
That is why the case r = is called the "no-go" one. It is characteristic for evolution equations 
in general. First the "no-go" case was completely investigated for the one-dimensional linear 
heat equation in [12]. It was proved that the problem of finding the conditional symmetry 
operators with the vanishing coefficient of dt is reduced to solving the initial equation. In [30] 
the proof was extended to the class of (1 + l)-dimensional evolution equations and in [22] this 
result was generalized for evolution equations with n space variables. 
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Let us note that "no-go" has to be treated as impossibility of exhaustive solving of the prob- 
lem. At the same time, imposing additional constraints on the coefficient 9, one can construct 
a number of particular examples of operators with r = and then apply them to finding exact 
solutions of the initial equation. It is the approach that was used in [15] for fast diffusion equa- 
tion ([3]). Since the determining equation has more independent variables and, therefore, more 
freedom degrees, it is more convenient often to guess a simple solution or a simple ansatz for 
the determining equation, which can give a parametric set of complicated solutions of the initial 
equation. (Similar situation is for Lie symmetries of first-order ordinary differential equations.) 

Consider the case r / which admits complete solving unlike the previous case. We can 
assume r = 1 up to the usual equivalence of reduction operators. Then the determining equations 
for the coefficients £ and 9 have the form 

( 8 ) 

9 X x — 99 x , &t — ^xv — ixx — £,xV + £,0 X — 99 v . 

Theorem 1. A complete list of Gi-inequivalent non-Lie reduction operators of the potential fast 
diffusion equation (|5j) is exhausted by the following ones: 

1. dt + sd x + f(co)d v , where uj = x + et; 

2- dt + f((jS){d x + d v ), where to = x + v; 

-2 

3. dt + Zdx + fot + cpxCjdv, where £= — ■ — , ip G {t + e x , tf(x)}; 

v + ip 

4. d t + td x -*i±*4d v , where £ = _ 2 1±2^, 

1 + x tan v — x 

X € {tan(2t) tanhx, coth(2t) cot x}; 

K o . ta Xt + Xxip. , t l-xtanhv 

5. d t + t,Ox : r^i where £ = -2- — , 

1 — x tann v — x 

( i / \ , . , , e 2:r tanh2t + l 2 - e 2x - e u ~\ 

X € <. tanh(2i) tanh x, tanh(2t) coth x, coth(2r) coth x, — =- , ^- jr \- 

L c — tann 2t 2 -\- & ~\~ e J 

Here e € {0, 1}, / is an arbitrary nonconstant solution of the ordinary differential equation 
fujuj = ffw) i- e - f £ { — 2/w, — 2cotw, — 2tanhw, —2 coth a;} mod G2. 

Proof. Here we only outline a sketch of proof. Any solution of the equation £„„ = belongs 
to the set {p, —2/(v + p), — 2/icotw, — 2/itanho;, — 2/icotho;}, where u> = fj,(v + p), [i and p> 
are arbitrary functions of t and x, /i 7^ 0. The second equation of dH) is a linear inhomogeneous 
second-order ordinary differential equation with respect to 9, where v is the independent vari- 
able and t and x are assumed parameters. It is possible to construct its partial exact solution 
without irrational singularities for any above value of £. The solutions of the corresponding 
homogeneous equation have irrational singularities if £ v 7^ 0. In view of the other equations 
of ©, the part of 9 containing such singularities has to vanish identically. Moreover, \i = const, 
i.e. n = 1 mod G2, and <p satisfies an overdetermined system of differential equation in t and x. 
Integration of it for all above values of £ and classification of obtained solutions up to equivalence 
with respect to G2 with excluding Lie cases result in the statement of the theorem. □ 

All operators from Theorem 1 are potential nonclassical symmetries of equation 
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5 Connection between classes of 

nonclassical and potential nonclassical symmetries 

Let us investigate connection between reduction operators of equations ([3|) and ([5]). 

As mentioned in the introduction, one of the problems studied in [15] was construction 
of partial classes of reduction operators for diffusion equations ([2]) in the "no-go" case when 
operators can be reduced to the form d x + rj(t, x, u)d u . Namely, Gandarias proposed to look for 
the coefficient rj with the ansatz 

rj l {t, x)u + rj 2 (t, x) 



/(«) 



(9) 



where f(u) = u~ a . After substituting ([9]) in the determining equation for r\ and splitting with 
respect to u, one obtains an overdetermined system for the functions tj 1 and rj 2 . For the fast 
diffusion equation ([3]) this system has the form 

222 22112,1 111 nn\ 

Vxx = V Vx, Vt = V Vx - V Vx + Vxx, Vt = V V x - ( 10 ) 



System ()10p can be derived from ([8]) with reduction by the group of translations with respect 
to v, i.e. with assuming that £ and 9 do not depend on v and re-denoting £ = —tj , 9 = rj 1 . 

This observation can be easily explained in a rigorous way for any pair of equations (pQ) 
and © with the same function /. 

Consider reduction operators Q = d% + £d x + 9d v and Q' = d x + r\d u of equations (|6|) and (JTJ) 
correspondingly, where the coefficients £ and 9 depend only on t and x, the coefficient r\ is defined 
by (J9|). The conditional invariance criterion applied to equation © (or ([T]) ) and the operator Q 
(Q 1 ) implies the following determining equation on £ and 6 (rj 1 and rj 2 ): 

(ttxVx 2 - (txO + co x )v x + ee x )f{v x ) 

+((6 + m x )v x -0 t - 29£ x )f(v x ) + {-ixxv x + O xx )(f(v x )) 2 = 

(or (r t 1 u + ri 2 )(rilu + r]l)f(u)- 

{(r,] - IrfriDu + rf t - 2r,° V 1 x )f(u) + (r^u + ^(/(u)) 2 = ) 

which has to be additionally split with respect to v x (u) . It is obvious that the systems obtained 
in the both cases after splitting coincide under the supposition r] 1 = rj 2 = 6. The charac- 
teristic equation Q[v] = 9 — vt — £v x = can be rewritten on the manifold of solutions of the 
potential system 

v x = u, v t = f(u)u x (11) 
in the form 

u x ji—— = 

fw 

and coincides in this way with the characteristic equation Q'[u] = 0. 
Therefore, the following proposition is true. 

Proposition 1. Q = dt + (,d x + 9d v , where £ = £(t,x) and 9 = 9(t,x), is a reduction operator 
of equation ([6]) iff 

Q =o x -\ d u 

fw 

is a reduction operator of equation (pTJ) . 

System (fTT|) establishes connection between the corresponding sets of invariant solutions. 
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6 Exact solutions of fast diffusion 
and nonlinear filtration equations 

All invariant solutions of ([3]) and ([5]), which were earlier constructed in closed forms with the 
classical Lie method, were collected e.g. in [21,25]. A complete list of Gi-inequivalent solutions 
of such type is exhausted by the following ones: 

1) u = — -r, v = - In \e~ x + ee*|; 

; l + se x+t 1 1 

2) u = e x , v = e x + t; 

1 , . . f d& 

3) u = — — — -77, v = ln\t\ + 



x-t + fite~ x / f ' " J i? - 1 + fie 



D=x/t 

X — t 



, 2t I It i X , 

4) u = — = v\ „ = , v\ ,=2 arctan — , v\ , = In 

' x 2 + et 2 ' le=0 x' le=1 t '^-i 

2t 

5) u = 2 — j v = 2itanx; 



x + 1 



cos^ X 
-2t 
cosh 2 a; ' 

2t 
sinh 2 x ' 



6) it = -g- — , u = — 2itanhx; 



7) it = - 2 — , v = — 2tcothx. 



Here e and are arbitrary constants, e € {—1,0,1} mod G\. The below arrows denote the 
possible transformations of solutions (|12p to each other by means of the potential hodograph 
transformation ([7]) up to translations with respect to x [25]: 

O 1) E =0 ; l)e=l < ► l)e=-l, a;+t<0 ! O l)e=-l, s+t>0 ! 2) < > 3)^=0, x>t ', 

O 4) e=0 ; 5) < — ► 4) e=4 ; 6) < — > 4) e= _ 4j ; 7) < — ► 4) e= _ 4j |^| >2 |* | . 

The sixth connection was known earlier [10,27]. If fi ^ solution 3) from list (|12j) is mapped 
by © to the solution 

8) it = - t + ixte~^\ uj = x - In 

which is invariant with respect to the algebra (tdt + d x + ud u ) . Here i? is the function determined 
implicitly by the formula f($ — 1 + [ie )~ 1 d$ = a;. 

Some classes of non-Lie exact solutions of ([3]) were obtained in [15,28,29]. These solutions 
and the ones similar to them can be represented uniformly over the complex field as compositions 
of two simple waves moving with the same "velocities" in opposite directions: 

a 2 

u = — (— cot(ax + (3t + 7) + cot(ox — fit + S)) 

(13) 

_ o2_ 2 sin(2ffl + 7 - 5) 

~~ cos(2/?t + 7 - 5) - cos(2ax + 7 + 5) ' 

where a, /3, 7 and <5 are complex constants, a/? 7^ 0. It can be proved that function (I13p takes 
real values (for real x and t) iff up to transformations from G\ 

(a, (3, 7, 5) £ {(1,1,0,0), (i,i,0,0), (M,tt/2,0), (t, i, tt/2, tt/2), (i, 1,0,0), (1,1,0,0)}. 
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Using representation (fl~3l) and the above values of tuples (a, f3, 7, 5), we obtain the following 
solutions of fast diffusion equation ^ and nonlinear filtration equation (0): 



1') it = cot (2; — t) — cot(x + t) = 
2') u = coth(x -t) - coth(x + t) = 

3') u = coth(x — t) — tanh(x + t) - 
4') u = tanh(x — t) — tanh(x + t) ■ 
5') u = cot (22 + t) — cot(ix — t) = 
6') u = i cot{x + it) — i cot(x — it) 



2 sin 2t 



v = In 



cos 2t — cos 2x ' 
2 sinh 2t 
cosh 2x — cosh 2t ' 

2 cosh 2t 
sinh 2x — sinh 2t ' 

2 sinh 2t 



sin(x — t) 



sin(x + 1) 

sinh(x — t) 



v = In 



u = In 



sinh(x + 1) 
sinh(x — t) 



cosh 2x + cosh 2t 
2 sin 2t 



u = In 



cosh(x + 1 
cosh(x — t) 



cosh 2x — cos 2t 
2 sinh 2t 
cosh 2t — cos 2x 



cosh(x + 1) 
, v = 2arctan(cott tanhx); 



, v = 2 arctan(cotht tanx). 



(All two's in the latter expressions for u can be moved over with scale transformations from G\.) 
Transformation (J7]) acts on the set of solutions l') _ 6') in the following way: 

l')cos2t<cos2x 4 * ^ / )\t^t+Tv/2,x^x/2,v^2v'i 1 )cos 2t>cos 2x < > 5')\x-w/2,v->2v— n i 



2' 



4') 



x|<|t| * 

O 3 / ) x< t| x ^ x /2, t) -^2t); 



/2,v^2v'i w 2 / )| a .| > | t || a ._ >a ./ 2 ,«-*2'y j 

3 ) x >t 



* 3')x>t\x-+-x/2,v^-2v : > O 6')\x-^x/2,v-^2v 



These actions can be interpreted in terms of actions of transformation ([7]) on the nonclassical 
symmetry operators which correspond to solutions l')-6'). 

In [29] Rosenau took advantage of additive separation of variables for the potential fast 
diffusion equation © and constructed solution 4'). Using the generalized conditional symmetry 
method, Qu [28] found solutions which can be written in forms 1') and 6'). After rectifying 
computations in two cases from [28], one can find also solutions 2') and 5'). Solutions 1'), 3') 
and 4') were obtained in [15], at least, in one from the above forms, but equivalence of these 
forms was not shown there. 

One of techniques which can be applied for finding the above solutions is reduction by con- 
ditional symmetry operators of the form Q = d x + (tj (t, x)u + r] 2 (t, x))ud u (see [15] for details). 
Namely, solutions 1'), 3') and 4') are obtained with the following operators: 

dx + (u 2 — 2 cot (a; — t)u)d u , d x + (u 2 — 2 coth(a; — t)u)d u and d x + (u 2 — 2 tanh(x — t)u)d u . 

We supplement the list of solutions adduced in [15, 28, 29] with similar ones, namely, with 2') 
and 5'). Solution 2') can be also constructed by reduction with the second above operator. Real 
solutions 5') and 6') correspond to the similar operators 

d x + (iu 2 — 2 coth(x — it)u)d u and d x — (iu 2 — 2i coth(t — ix)u)d u 

with complex- valued coefficients. 

All reductions performed with reduction operators from Theorem 1 or with equivalent ones 
result in solutions which are equivalent to the listed Lie solutions 1)— 7) or solutions l')-6'). For 
example, the operators 

dt — d x — 2 cot (a; — t)d v , dt — d x — 2 coth(a; — t)d v and dt — d x — 2 tanh(x — t)d v 

lead to solutions 1'), 3') and 4') correspondingly (see also Section [S]). 
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7 Conclusion 



In this Letter we present classification of reduction operators for nonlinear filtration equation §5§ 
with summary of necessary notions and statements, a basic sketch of the proof and a list of 
constructed exact solutions including both Lie and non-Lie ones. Since (0) is the potential 
equation of fast diffusion equation ([3]), all the obtained operators are potential nonclassical 
symmetries of (|3j). Moreover, most of them are nonprojectible on the space of the independent 
variables t and x that leads to technically cumbersome implicit reductions of © to ordinary 
differential equations. Now we optimize the proof and hope to realize it in a quite compact form. 
Presentation of the proof and reduction technics will be subjects of a future paper. 

We continue our investigation on potential reduction operators of the nonlinear diffusion 
equations from class (pQ). In some sense, equation ([3]) is singular in this class with the potential 
nonclassical symmetry point of view. More precisely, as a result of joint work with Prof. Sopho- 
cleous the following theorem have been proved recently. 

Theorem 2. Nonlinear filtration equations © admit non-Lie reduction operators with non- 
vanishing coefficients of dt only in the case of the Fujita's nonlinearities 



where a, b and c are constants. 

Let us note that there are exactly three G~-inequivalent cases of the Fujita's nonlinearities: 



The equivalence group G~ of class ([6]) is formed by the transformations 

t = £l t + e 2 , x = e[x + e' 2 v + e' 3 , v = e'{x + e'^v + e 3 ', f = ej 1 (e' x + e' 2 v x ) 2 f , 

where ex, £2, e^, e" (i = 1,2,3) are arbitrary constants, £x(£x e 2 ~~ £ 2 £ i) 0- The nonclassical 
(conditional) symmetries of the (l + l)-dimensional linear heat equation (/ = 1) were completely 
studied in [12]. Analogous investigation of the second case (/ = v" 1 ) is carried out in this Letter. 
Therefore, to complete classification of reduction operators in the class of nonlinear filtration 
equations ([6]) with respect to G~ (see Definition [3]) , it is enough to describe reduction operators 
of the equation with the latter nonlinearity, and we achieved significant progress in solving this 
problem. 
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